Classical projected phase space density of billiards and its relation to the quantum 

Neumann spectrum 
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A comparison of classical and quantum evolution usually involves a quasi-probability distribution 
as a quantum analogue of the classical phase space distribution. In an alternate approach that 
we adopt here, the classical density is projected on to the configuration space. We show that for 
billiards, the eigenfunctions of the coarse-grained projected classical evolution operator are identical 
to a first approximation to the quantum Neumann eigenfunctions. However, even though there 
exists a correspondence between the respective eigenvalues, their time evolutions differ. This is 
demonstrated numerically for the stadium and lemon shaped billiards. 



A comparison of classical and quantum dynamics in 
terms of appropriate evolution operators is of much in- 
terest [1-4] . Besides improving our understanding of the 
statistical properties of quantum spectral fluctuations, 
these studies are expected to shed light on the effective 
irreversibility (e.g. relaxation to the invariant density) 
observed in low-dimensional isolated systems which are 
otherwise time reversible. 

In the probabilistic approach to dynamics, classical 
time evolution is studied in terms of the propagation of 
the phase space density, p(p,q). The evolution of p(q,p) 
is governed by the Perron- Frobenius (PF) operator, 



£} o p(x) = / 5(x - f*(x'))p(x')dx' 



(1) 



where x = (q, p) is a point in phase space and f*(x) is 
its position at time t. In the Hilbert space of phase space 
functions, £* is unitary [5,6]. Thus, the eigenvalues lie 
on the unit circle. The existence of an invariant density, 
poj implies £* o p = p . A system is ergodic if the unit 
eigenvalue is nondegenerate. A knowledge of the spectral 
decomposition of £* allows one to evaluate correlations, 
averages and other quantities of interest. For integrable 
systems, the spectrum of £* is discrete while for mixing 
systems, there is a continuous spectrum apart from the 
unit eigenvalue. The decay to the invariant density is 
connected to the continuous part of the spectrum. When 
the decay is exponential as in case of a class of chaotic 
systems, the fourier transforms of time-correlations have 
poles in the complex frequency plane or broad peaks (res- 
onances) along the real frequency axis, the positions of 
which arc independent of the observable chosen [7] . 

A starting point for the comparison of classical and 
quantum dynamics is usually a quasi-probability distri- 
bution involving the density operator, p = \ip >< ip\. 
This enables one to "lift" the quantum state to the phase 
space. The process is however not unique and depends on 
the ordering scheme used to arrange the non-commuting 
operators (q,p) [8,9]. A commonly used quasi-probability 
distribution is the Husimi function which is a coherent 
(most classical) state representation of the quantum den- 
sity operator. The Husimi propagator is thus a quantum 
analogue of the Perron- Frobenius operator and the eigen- 



states of the two provide a means of comparing quantum 
and classical evolution. 

While the exact Husimi function does not decay to the 
invariant density, it is found that coarse graining of the 
phase space leads to a loss of unitarity [3,4]. For the 
kicked top [3] and maps on the torus [4] , the eigenvalues 
of the coarse grained classical and quantum propagator 
have been found to be identical. 

The primary aim of this paper is to explore the ex- 
istence of a correspondence between the classical and 
quantum spectrum for the class of systems referred to 
as billiards. However, rather than "lifting" the quantum 
description to the phase space using a quasi-probability 
distribution, we shall project the classical density to the 
configuration space by integrating out the momentum: 
p(q) = J p(q, p) dp. While we shall restrict ourselves to 
billiards here, the essential idea of seeking a quantum- 
classical correspondence using the projected density in- 
stead of the full phase space density is applicable to other 
systems as well. In the following, we shall show that 
for billiards, there exists a correspondence between the 
eigenvalues of the coarse-grained projected classical evo- 
lution operator (£p) and the quantum Neumann spec- 
trum, while the respective eigenfunctions are identical to 
a first approximation. However, the eigenvalues evolve 
differently with time so that classical and quantum evo- 
lution differ [10] 

Apart from being a paradigm in the field of classi- 
cal and quantum chaos, billiards have relevance in vari- 
ety of contexts. The Helmholtz equation describing the 
quantum billiard problem also describes acoustic waves, 
modes in microwave cavities and has relevance in stud- 
ies on "quantum wells" , "quantum corrals" , mesoscopic 
systems and nanostructured materials. 

In a classical billiard, a particle moves freely inside a 
given enclosure and reflects specularly from the wall at 
the boundary. Depending on the shape of the bound- 
ary, billiards exhibit the entire range of behaviour ob- 
served in other dynamical systems. They also provide a 
means of coarse graining that is perhaps unique. The 
dynamics of billiards with smooth boundaries can be 
coarse grained by polygonalizing the boundary [13,14,12] . 
Rational polygonal billiards are non-ergodic and non- 
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mixing. However, the short time dynamics of a polyg- 
onalizcd billiard can approximate that of the smooth bil- 
liard [13]. 

The quantum billiard problem consists of determin- 
ing the eigenvalues and eigenfunctions of the Hclmholtz 
equation V 2 ip(q) + k 2 ip{q) = with ip(q) = (Dirich- 
let boundary condition) or n.V^ = (Neumann bound- 
ary condition; n is the unit normal) on the billiard 
boundary. Its semiclassical description holds the key to 
the quantum-classical correspondence. However, since 
a quantum state (or the quasiprobability distributions 
constructed out of it) can essentially resolve phase space 
structures of the size of a Plank cell, polygonalization 
provides just as much information about the quantum 
state at the semiclassical level [14,15]. 

We shall consider a polygonalized billiard as an ex- 
ample of a coarse-grained system. The "unfolded" dy- 
namics of a polygonalized billiard can be viewed locally 
as a straight line on a singly connected invariant surface 
consisting of multiple copies of the enclosure glued to- 
gether appropriately at the edges, each copy denoting a 
momentum direction that is related to the previous one 
by the law of reflection at the glued edge (for a figure, 
see [11,19]). As the magnitude of the momentum (p) and 
the angle <p that p makes with (say) the X-axis are con- 
served, it is convenient to treat p in polar coordinates 
(p, ip). Transforming from (p x ,Py) to (p, ip), the p inte- 
gration in cq. (1) simplifies as 

J dp x dp y 5(p x -p' x t (q',p')S(p y -p'y{q[,p')h(q!,p') 
= J dip5(ip-v')h(q! u ,<p'-,p) = h(q[ u ,<p;p) (2) 

where q u is a point on the unfolded space, p^,'(q',p') 
(p' y ) is the x (y) component of the momentum at time t 
for the initial phase space coordinate (q',p') and 

K<h,wp) = ^(qu - qu"(qu';¥',p))p(qu;¥',p) (3) 

Thus 

£ t op = J dq' u <5(q u - q u '*(q' u ; ip,p))p(q[ u ). (4) 

Note that the time evolution of p depends on p through 
the kernel as q u " depends on both the initial position 
and momentum. Thus £* = C l {p). 

Projection onto to the configuration space requires an 
integration over the angle p>. The time evolution of the 
projected density is thus given by: 

1 f 2n 

4°p(q) = ^ y o dpC t (p)op( q ) 

= 2^ / dc[u ' dip 5 ( <lu -qu'*(qu';p>^)Mq') ( 5 ) 

The spectrum of £ p , can be studied by evaluating its 
trace 



Tr£* P = i-/#]T e ^ (6) 

= 2~7r J dqu f dip S ( q - u ~ qu'*(qu';p,y))- (7) 

Note that due to the multiplicative nature of the Perron- 
Frobenius operator, its eigenvalues A n (i; <p) are of the 
form e A "(^*. It can be shown that for t > [16,11] 

OO 

Tr L t P ~NC + NY J 9{VE~nl) (8) 

n=l 

where {E n } refers to the quantum Neumann spec- 
trum, I = tv where v refers to the velocity, g(x) = 
(nx) cos(x — 7r/4), N is the maximum number of 
allowed momentum directions and C ~ 1/N is a con- 
stant [17,18]. Since g(x) ~ ^ e ix sin ^ dip for large 
x, it follows that for v = 1, 

\ n (ip) = i^/E^am(p) (9) 

Thus, the power spectrum of a projected density contains 
peaks at \J E n . 

The above correspondence between the classical and 
quantum spectrum of polygonal billiards arises from a 
similarity in the traces of C} p and the quantum propaga- 
tor when expressed in terms of periodic orbits. The cor- 
respondence however strictly holds for large E n so long 
as one uses a delta function kernel in £* [11]. For smaller 
value of E n , the correspondence exists if the delta func- 
tion kernel in cq. (5) is smoothened [19,11]. This effec- 
tively results in a coarse graining of the dynamics and in 
generic cases, leads to the the inclusion of a higher order 
term in the classical and quantum trace [11,19]. Despite 
the coarse graining, {E n } in most cases refers to the 
approximate quantum Neumann spectrum as the semi- 
classical trace formula, which is used in arriving at the 
correspondence, remains inexact. For integrable poly- 
gons such as the rectangle, the correspondence between 
the spectrum of the projected Perron-Frobenius operator 
and the Neumann spectrum is exact and can be shown 
directly [11]. 

Note that despite the correspondence, quantum time 
evolution differs from evolution due to C l P as the eigen- 
values evolve differently The quantum Neumann eigen- 
functions are however approximate eigenfunctions of C P . 
This has been established [12] for a quasiclassical adap- 
tation of C P when the Dirichlet eigenstates are of in- 
terest. For the Neumann problem, a similar derivation 
follows provided the quasiclassical kernel is replaced by 
the classical kernel. We shall demonstrate numerically 
that V'n(q) indeed approximates the quantum Neumann 
eigenfunctions. 

For generic polygons, the unit eigenvalue of £* is non- 
degenerate and the corresponding eigenfunction is a con- 
stant. This is true as well for the projected operator, C p . 
Thus Ao = 0. The quantum Neumann problem also has 
a constant as its ground state eigenfunction and the cor- 
responding eigenenergy Eq = 0. This is consistent with 
the results presented here. 
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In order to determine the eigenvalues and eigenfunc- 
tions of £p, we shall first evaluate its smoothened kernel 

A- P (q,q',t) = i- J #^(q-q'*M) 

= ^Vn(q)V-n*(q')An(i) (10) 



as a function of time. Here 5 e is a smoothened delta func- 
tion and V'n(q) are the eigenfunctions of the projected 
Perron- Frobenius operator, £ P . As an example of the 
smoothened delta function, we consider the hat function 
which is zero outside a cell of size e [20] . The ip integra- 
tion is performed by shooting trajectories from a point q' 
at various angles and evaluating the fraction of trajecto- 
ries in a cell of size e at q [11]. Since A„ = i^/E^ sa\{Lp) , 
for v = 1, a fourier transform of Kp(q, q', t) has peaks 
at k = \fEm, the width depending on e and the heights 
on Vrt(q)- An eigenfunction corresponding to a partic- 
ular eigenvalue can thus be measured by varying q and 
measuring the height of the desired peak. 
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FIG. 1. An eigenfunction of Cp corresponding to the first 
non-zero eigenvalue of the polygonalized stadium. It is sym- 
metric in Y (Neumann) and antisymmetric in X (Dirichlet). 
In the region outside the stadium, the value has been set to 
zero, (b) its quantum counterpart in the smooth stadium. 

We now present our numerical results for the stadium 
and lemon shaped billiards. The chaotic stadium shaped 
billiard that we consider, consists of two parallel straight 
segments of length 2 joined on either end by a semicir- 
cle of unit radius. This has been polygonalized using 



10 segments to approximate the semicircle. Figure la 
shows the eigenfunction of C P corresponding to the first 
peak at a non-zero k (k ~ 0.70) in the power spectrum 
of (the smoothened) kernel Kp(q, q, t) for the polygo- 
nalized stadium. Only the first quadrant is shown due 
to the reflection symmetry of the system. Note that we 
have plotted the intensities as the peak heights are pro- 
portional to \ip n (q)\ 2 - Figure lb shows the correspond- 
ing quantum Neumann eigenfunction of the smooth sta- 
dium at k ~ 0.87 found using the boundary integral 
technique. The eigenfunction of the projected Perron- 
Frobenius operator clearly approximates the quantum 
Neumann eigenfunction. This is true for other stadium 
eigenfunctions as well. 
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FIG. 2. (a) An eigenfunction of C P corresponding to the 
third non-zero eigenvalue of the polygonalized lemon shaped 
billiard. It is symmetric in X and antisymmetric in Y. (b) 
its quantum counterpart in the smooth lemon with k ~ 3.51. 

We consider next a lemon shaped enclosure con- 
structed by the intersection of two circles of radius 2.5 
centred at (2,0) and (-2,0) respectively. Each arc is ap- 
proximated by seven segments. Figure 2 shows a compar- 
ison similar to fig. 1 for this enclosure for the third peak 
in the power spectrum of Kp at k ~ 3.08. The agreement 
with the exact quantum eigenfunction (k ~ 3.51) of the 
smooth lemon shaped billiard is again reasonable. 

We have thus seen that there is a correspondence, al- 
beit approximate, between the eigenstates of the quan- 
tum and projected classical evolution operators. The 
theoretical basis (see also [16,11]) clearly indicates that 
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the eigenstates obtained using this method are at best 
"semiclassical" in nature since the quantum trace for- 
mula, which is necessary in order to connect the eigen- 
values of Hp with the quantum Neumann eigenvalues, is 
only approximate as higher order corrections have been 
neglected. In cases where the corrections are zero (such 
as in the rectangular or equilateral billiards), the quan- 
tum Neumann eigenstates are identical to the eigenstates 
of Cp and the smoothcning parameter e can be made 
small. In general, at higher energies, the peaks in the 
fourier transform of Kp(q, q',t) are generally harder to 
resolve as the density of eigenvalues increases with k. 
Thus, measuring projected cigenfunctions using trajec- 
tories becomes harder. 
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proach to the uniform density, p av , is highly oscillatory 
and the overall decay rate may differ from t -1 / 2 due to 
the sum of oscillatory terms. We have estimated the over- 
all decay by locally averaging p(q, t) over time. For each 
of the three cases in fig. 3, the best fit of a + b/t 1 / 2 is also 
shown. These results testify that there is a decay to the 
uniform density in polygonalized enclosures. 

Finally, a discussion on coarse-graining is important 
to understand the significance of the result. There are 
two levels at which this has been carried out. First, the 
boundary has been polygonalized. This enables us to 
connect the quantum and classical eigenvalues analyti- 
cally. As the number of segments can be increased to ap- 
proximate the smooth billiard arbitrarily well, it might 
be expected that the eigenvalues and eigenfunctions con- 
verge to those of the smooth billiard in both the quantum 
and classical case. We have verified this numerically for 
a few eigenstates of stadium and lemon billiards. Thus 
coarse graining via polygonalization might be unimpor- 
tant for observing the correspondence. The second and 
more significant coarse graining of the dynamics is related 
to the smoothcning of the delta function kernel. This en- 
ables us to connect the classical and quantum eigenvalues 
even for smaller values of A„ and is hence indispensable. 

Finally, it is worth noting that the determination of 
exact eigenstates of the Perron- Frobenius operator is gen- 
erally nontrivial while quantum states are easier to deter- 
mine. The usefulness of the results presented here thus 
lies in using quantum states to evolve classical configu- 
ration space densities. 



FIG. 3. The time evolution of the locally space and time av- 
eraged configuration space density, p(q,t), in a lemon shaped 
billiard with different degrees of polygonalization. The num- 
ber of linear segments approximating each arc is shown in the 
figure. The point q is different in each case. The background 
shows the density at t ~ 1250 for the 31 segment case. 

Note that there are basic differences in the interpreta- 
tions of the quantum and projected classical eigenstates. 
Each quantum eigenfunction is associated with a density 
which is invariant in time while only the constant eigen- 
function of C P qualifies as a configuration space density 
and is invariant. All other densities, p(q), decay to the 
invariant density on evolving with C p . This is true even 
when the dynamics is integrable. The eigenvalues of C P 
thus form a decay spectrum. As an example of the de- 
cay to the uniform density in polygonalized billiards, we 
present in fig. 3, the evolution of a localized projected 
phase space density, locally averaged in space and time, 
for three different versions of the polygonalized lemon- 
shaped billiard considered above. The time evolution 
of a projected density can be expressed in terms of the 
eigenstates of C P as 

p(q, t) ~ p av + c nV"n(q) COs(V£^ - 7T/4) (11) 

n=l 

for t sufficiently large. The coefficients, {c„}, apart from 
some factors, depends on J ^> n (q)p(q, 0)dq. The ap- 
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